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MASS SPECTROMETRIC. ANALYSES OF HYDROCARBON. 
AND GAS .MIXTURES' 


By A. Keith Brewer and Vernon H. Dibeler 


ABSTRACT 
The basic. prmeiples underlying mass spectrometric analysis. of hydrocarbon, 
and gas mixtures are outlined. The method of calculating the composition from 
e mass-abundance records is described. Analyses of 4 number of hydrocarbon 
nd other gas mixtures are shown. Whenever possible, comparisons with different 
hods of analyses are given. The reproducibility and the accuracy obtainable 
Ass spectrometric analyses are described. 
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I. INTRODUCTION 


The basic principles involved in the mass spectrometer were worked 

t as early as 1910 by Sir J. J. Thomson. The first mass spectro- 
rraphs utilizing these principles, much as in ‘the present-day instru- 
meat, were built by Aston and by Dempster about 1920. In spite 
of this early conception, it is only within the last 2 years that the mass 
spectrometer has been developed .to the point where it is capable of 

aking precise analyses of gas: and hydrocarbon mixtures. 

T his new development has been made possible, not from changes in 
principle or basic design, but from the discovery thiat hydroc arbon 
nol cules ae bombardment at low pressures (107° to 10-* mm) by 
electrons having 50 or more volts of energy are dissociated into all 
possible disintegration fragments.?. The ratios in which the fragments 
occur are the same over a wide pressure range for ‘each molecular 
species but are never the same for different species. The development 
of sensitive electronic pick-ups and automatic-recording mechanisms 
has not only contributed to the ease of operation. but has made possible 
a speed and precision that could not be attained by manual controls. 


rk was financed in part from funds made available by the Rubber Reserve Co 
ver, Jr. and H. W. Washburn, Calif. Oil World Petroletim Ind., (Nov. 1941); HW. W. Washburn, 
] und S. M. Rock, Ind. Eng. Chem. Anal. Ed. 15, 5 1943); H. W. Washburn, H. F. Wiley, 


«' he 
ck, and C. EF. Berry, Ind. Eng. Chem. Anal. Ed. 17, 74 (195 
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This paper describes briefly the method by which chemical com. 
position can be determined from mass spectrometric records, the types 
of mixtures that are being analyzed;.and the reproducibility that js 
being realized. Representative analyses made in the mass spectrom. 
eter laboratory at the National Bureau of Standards by the use of , 
mass spectrometer manufactured by the Consolidated Engineering 
Corporation are cited to illustrate the applicability of the instrument, 


II. DESCRIPTION OF THE MASS SPECTROMETER 


In the mass spectrometer, molecular masses are separated by virtue 
of differences in their momenta when all have been given an identica| 
kinetic energy in a given direction: Four steps are involved in the 
separation. 

First, the molecules are given a positive charge. This is accon.- 
plished in the instrument by bombardment with a stream of 50-voli 
electrons, the electron current being about 5 microamperes. The gas 
pressure in the ionizing chamber is of the order of 10-° millimeter 
The filament is placed outside the ionization chamber, where the ga 
pressure,.as measured on an ionization gage, ranges from 107? to 1)" 
millimeter. . This insures long filament life and prevents products of 
molecules cracked on the hot filament from reentering the ioniziag 
chamber. 

Second, the ions are given an identical kinetic energy along a given 
path by a potential difference between two accelerating slits. In the 
instrument this potential can be manually or automatically varied 
from 4,000 to 200 volts.. The kinetic energy of the ions is given by 
1/2 mv’=eE, where m is the ion mass, v the velocity, e the charge on 
the ion, and E the accelerating potential. 

Third, the ions are resolved in a uniform magnetic filed by being 
deflected in arcs whose radii are proportional to their momenta. - Thy 
instrument is of the 180-degree focusing type, the entire resolved ion 
path being within the uniform magnetic field of a large electro magnet 

The centripetal force deflecting the ions is balanced by centrifugal 
force, as expressed by Hev=mv’/r, where H is magnetic field strength 
and ris the radius of deflection. Dividing both sides by »v and re- 
arranging, the above expression gives my= Her; this shows the radius 
to be proportional to the momentum. Substituting the energy valu 
for v in the momentum equation, the basic expression for the principl: 
of the mass spectrometer is obtained ‘as (m/e) =(H?r*)/2E. Express 
ing all values in practical units this equation becomes 


«3 10"=144.00¥m7/V, 


21.66 10-* xk mV}, 
Hr: [? 56 m 


4.803 X10-" 300, 


where H/ is in gauss, V in volts, r in centimeters, end m in units of 
hydrogen atoms. 

Fourth, the resolved ion beams are collected and recorded. In the 
instrument the focusing radius, r, is fixed at. 12.7 centimeters. The 
various masses are made to pass over the collector slit by varying | 
while H is kept constant. This is accomplished by permitting the 
charge on the ion-accelerating slit to leak to ground through a high 
resistance: As the potential decreases 9 progression of masses fro 
light to heavy is brought to. focus on the collector slit. The arrange 
ment chosen has the advantage of providing a nearly linear record 8 
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the radius of the ion path varies as the{square root of the mass, 
whereas the rate of potential leak is inversely proportional to the 
voltage. 

The resolved ion masses dre collected and the ion current is amplified 
by a direct-current feed-back amplifier. The output is fed through 
four galvanometers ‘in series. Shunts are provided to reduce the 
sensitivity of the second, third, and fourth. galvanometers to one- 
third, one-tenth, and one-thirtieth that of the first. This gives the 
g-inch-wide recording paper an effective width of 240 inches. The 
deflections of the four galvanometers as well as thie scale divisions are 
photographed simultaneously on the moving record. , 

A mass marker is provided to make possible the easy identification 
of the various mass numbers. The marker consists of a relay ‘that 
trips at definite intervals to throw a‘spot of light on the moving 
record when the acceleration potential passes fixed multiples of a 


‘chosen reference voltage. 


III. MASS SPECTRA.OF HYDROCARBON MOLECULES 


The ionization potentials of hydrocarbon molecules range from 
about 9 to 14 volts. The 50-volt ionizing electrons possess enough 
energy not only to ionize the molecules but also to dissociate them 
into all possible fragments. As sufficient molecules are ionized for the 
probability laws to hold, the fragments of each molecular species 
always occur in the same abundance ratio. 

In-table 1 are given the cracking patterns for some of the-C, hydro- 
carbons. ) 

TaBLe 1.—Mass spectrometer calibration patterns 





| | 


n-Butane Isobutane} Butene-1 





| 1,2-Bu- | 1,3-Bu- | 
| tadiene | tadiene 


| 
cis-Bu- | trane-Bu- — 
|*s J 





0. - 0.44 
2. 4: Lan 

14. ; 24. 70 

40. 33. 70 

j . 50 


658. £ 
99. 
1444 
1208 
3666 
118.7 
1.41 | 


Trace 


23. 03 

5. 45 

20.00 

3. 23 
16. 57. | 

13. 04 

110. 30 

| 100. 00 
35 4.35 
. 07 .07 | 


a 
KP CMO, 











14. 01 19. 11 | 16. 76 
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In preparing table 1 the abundance of the base mass was assign 
an arbitrary value of 100; and the abundance of all other mass 


weighted accordingly. The base mass is that of the normal mole 


composed of C™ and H atoms. The only heavier peaks observed ¢ 


those due to the heavy isotopes, C™ and D.. No association betwe 


molecules -is detected at the pressure existing in the analy: 
chamber. 
The sensitivity for each compound is given at the bottem oi 


table. Sensitivity is defined as the number of divisions of base pe: 


per micron of gas pressure. As the peaks are proportional to 
pressure over the range in- which the instrument is operate: 


fensitivity is independent of. the sample pressure or of the part 


pressure of any constituent present. 
The isotopic contribution to successive peaks is shown in.-table ‘ 


TABLE 2 I soto pe correction factors 


“Hs . oth 
sHe 57 
:H ‘ 58 


4 
4 
| 
‘ 
4 
1 
+ 
t 
4 
4 
1 


The.C® and D contributions for the ions ‘given ‘in column 1 (tabi 


2) to the peaks 1 and 2 mass units higher (m/e+-1 and m/e 


computed on the basis of the observed abundance of C™ as 1.05 


percent and D as 0.02 percent. 
The masses given in Table’l-are due to singly charged ions 
addition to. these, peaks are observed for doubly. charged ions 


ghosts, and fer “labile” groups, principally hydrogen. In the: cas 


of n-butane, for instance, the singly charged ions constitute 


percent of the total. The percentage of doubly charged ions is highe 


for unsaturated hydrocarbons and: for diatomic molecules. 

Only singly charged ions are used. in eomputing the composit 
from a mass record. The presence of other ion groups introduces 
complication as they occur.in proportion to the pressure and are 
distinguished. They are in-fact-often a help in hedalaieg 
parent molecule. 

A mass record for a typical gas mixture is shown in figure 1. 
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The record has-been illuminated to show the contributions to th 
various mass numbers made by each constituent. The traces of th 
four. galvanometers decrease in sensitivity from top to bottom. 1] 
obtain the abundance of the various masses it is necessary to correc 
the observed peak height for the galvanometer tracing used. Tly 
highest galvanometer peak on thé paper is always selected, and th 
number of divisions it. subtends is multiplied by the shunt in order ty 
express the peak heights in terms of the most sensitive galvanometer 


METHOD OF CALCULATION 


Both the base peaks and the ratios of the various fragments in the 
cracking patterns must be considered in analyzing 9 mass spectrogran 
like that shown in figure 1. In consequence, it is necessary to have the 
cracking pattern and the sensitivity of each compound that might bk 
present before making the calculation. 

The procedure can best be illustrated by considering a simple hydro. 
carbon mixture having the pattern given in column 2 of table 3. The 
values given represent the relative abundance of the various ions 
recorded at each mass listed in column 1.. The records are all read 
0.1 division. All masses observed are recorded, except those ‘beiow 
24, which are not essential to the calculation when methane, nitrogen 
and carbon monoxide are known to be absent. 


TABLE 3.—Mézture analysis 








| 


| | 
Mixture |n-Butane Isobutane! 1-Butene | 2-Butene om, 


1, 3-Bu- | 1, 2-Bu- | 


tadiene | tadiene | 





0. 05 | 
. 52 
6. 00 

18.26 
3. 14 
. 84 
& 





-ll 


Caw wNoODa~ 





ofr @ 

















Ne oDoos@NQ 








| 


* S=sensitivity: «micron of pressure. 


Several methods can be used in separating the various components 
of the gas mixture. The one described herein is the principal proce- 
dure used at the National Bureau of Standards. It consists in re 
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moving mathematically the contribution of the various constituents 
from each mass number. When possible each compound is subtracted 
individually starting with the highést masses. ; 

A survey of the pattern in column 2 shows clearly that no hydro- 
carbons are present heavier than butanes of base mass 58; the 59 peak 
can be accounted for fully. by the heavy isotope contribution from 58. 
The first problem is, therefore, to determine if 58 is due to n-butane, 
to isobutane, or to both. Table 1 shows that n- and isobutanes make 
very different contributions to the 43 peak. As the 44 peak is the 
correct size for the isotope contribution from 43, propane and carbon 
dioxide are shown to be absent. In consequence, the only contributors 
to 58 and 43, other than butanes, ‘are the second isotopes from 56 and 
41; these contributions are 0.05. and 0.04 division, respectively. 

The relative contribution of n- and isobutanes to the 58 and 43 
peaks are expressed by the simple linear equations as- : 


43 peak=7.3476 n-butane+36.6667 isobutane=62.35 divisio 


58 peak=1,.00 n-butane+ 1.00 isobutane = 3.36 divisions. 

These can be solved by determinants as follows: 
162.35 36.67 | 
5. SE”. a 08 divisions 

7.347 36.67 ~—29.3191 ~~ : 

1.00 1.00 | 
Similarly, isobutane= 1.28 divisions. . These values appear in columns 
3and 4 at mass 58. When the two base peaks are multiplied by their 
respective patterns the contribution of the two butanes to all the 
masses is obtained and can be subtracted out. 

The second step is to remove the butanes shown to be present by 
the large 56 peak. For sake of simplicity 2-butene will be considered 
as an equimolar mixture of cis- and trans-isomers. The procedure is 
the same as that followed for the. butanes, but three equations are 
needed as three butenes are possible. Table 1 shows that with .refer- 
ence to the base peak, the peaks at 55 and 41 are markedly different 
for the butenes. The only possible remaining compounds tha can 
contribute to the 55 and_56‘peaks are the isotopes from 1, 3-butadiene 
and ethyl acetylene. Possible remaining contributors: to 41 are 
propenes and the isotope from methyl] acetylene. 

The approximate nonbutene contributor{to'the three peaks is first 
determined; this is later. corrected to fit the equations. Propene, 
having a base peak of 42, can be determined quite accurately on the -- 
first approximation. The butene contribution to 42 is small and does 
not differ greatly for the three butenes. Taking an average contribu- 
tion for the butenes, the size of the 42 peak ‘can be calculated from'the 
56 peak. In the present case the entire 42 peak can be accounted ‘for 
by butenes, hence only a trace of propene"can be’present. » The*fact 
that the 39 peak is almost a8 large as‘54 indicates ‘a*large concentra- 
tion of 1, 3-butadiene and'a negligible concentration‘of ethyl acetylene. 
The isotope contributions to the.41, 55, and 56 peaks can now be 
computed from table 1; the values are 1.37, 45.71 and 0.74 divisions, 
respectively. . 


n-Butane= 


654385452 
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The butenes are now separated by the following linear equations 


2.2819B, +1.6411B,+1.8723B,=117.86 divisions. 
0.4069B, +0.3754B,+-0.3128B,= 23.83 divisions. 
1.00B,. +1.00B, 1.00B, = 61.75 divisions. 


41 peak 
55 peak 
56 peak 


Solving by determinants; 


Butene-1= 
1117.86 
23.83 
61.75 
2.2819 
0.4069 
1.00 


1.6411 i. 
0.3754 0.2 
1.00 1.18 


0.0472 


1.00 


1.6411 


0.3754 0 
1.00 l 


8723 
3128 
00 


8723 
3128 


25.0 divisions. 


Similarly, butene-2 =34.5 division and isobutene=2.3 divisions. 

The above values for the butenes appear in columns 5, 6, and 7 
mass 56. The contribution of each butene can now be computed a 
subtracted. The remaining divisions are due to the butadienes a 
to possible smali amounts of methyl or ethyl or vinyl acetylene a 
to propene. Ethyl acetylene can be shown to be present jn no nm 
than a trace as the 38 peak is 7.84 percent of the 54 peak; this 
almost exactly the pattern for 1,3- butadiene. The ethyl acetylur 
pattern is almost twice this size. 

The butadienes can now be resolved in the manner’ described | 
the butanes and the butenes by using the 54 base peak and the 39 pea 

Thus 


0.9525Bdj,5- +0.3915Bdi,2 
L.00Bd,,; +-1 OOBd,,» 


1,3-Butadiene= 1053.65 divisions 
1,2-Butadiene=5.25 divisions. 


- 1005.66 divisions 
1058.90 divisions 


39 peak= 
54 peak 


These values appear.in columns 8 and 9 (table 3) at’: mass 54. 
contributions of the butadrenes can now be computed and subtract 
The residuals are given-in eolumn 10. 

The residuals are particularly important in that they show: 
in approximatrons as well as errors in readings and in computations 
Positive residuals indicate. too little estimated, whereas nega‘ 
residuals indicate. too much subtracted: . In the. present 
residuals ate all-small and within the limits of reading error. 
indicates .the correctness of the assumption that propene an 
various acetylenes were absent. 

The accuracy of the above assumptions can be estimated ~ 
sidering the size of the residuals that would have been left had 1 p 
cent of the various possible contaminants been present. These va 
would be: for propene, 18 divisions on 42; for carbon dioxide, 20 
sions on 44; for ethyl acetylene, 19 divisions on 54; for vinyl acetyle! 
21 divisions on 52; and for propadiene, 11 divisions on 40. In ad 
tion, each compound would contribute to the residuals on peaks 
lower mass than those mentioned. The residuals on all peaks 
weighed in the final calculations. 

In order to calculate the mole. percent of each constituent, the nu 
ber of divisions of the various components must be divided by | 
sensitivity. ‘This converts divisions to microns of pressure. The sw 


( 
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of all the partial pressures gives the total ‘‘computed”’. pressure of the 
sample. This pressure should be identical to the “read” pressure as 
riven by the mercury manometer on the inlet measuring volume. The 
read and computed pressures, however, usually differ by a few tenths 


Ci 
illimeter, due primarily to water vapor in the sample, which is 


a@ mili 
sorbed on and later desorbed from the dry glass walls of the inlet 
n. As water is almost always present and cannot be measured 
weurately, it is never advisable to compute the concentration of any. 
omponent by difference. All percentages, therefore, are figured by 
the use of the ‘“camputed”’ microns as a measure of total pressure of 

sample. , 5 


V. REPRESENTATIVE ANALYSES 


The types of gas mixtures that can be analyzed with the mass 
spectrometer can best be illustrated by representative analyses made 
t the. Bureau. 

The results obtained from a series of butadiene- samples are shown 

table 4. Included in the table are the freezing point values for 

-butadiene as measured by I. D. Rossini. 


TABLE 4. Vass Spectrometer (n ole percent)*® 


le 6 : pl Sample 8 ample { ‘Sample 10 


,-+-(), 04 99. 474.0. 03 97. 974 98. 57+0. 01 
O4 09+. 03 7 5 


FREEZING POINT (MOLE PERCENT) > 


Os t 99.16 +0.08 | 99.40 +0. 08 7.99 +0 22 +0. 08 | 98.10 


esent maximum deviation from average of 2 


ertainty estimated by F. D. Rossini 
The reproducibility of independent, analyses is illustrated in table 5 
TaBLe 5.~ Refined butadienes 


Cylinder 11-20 


1141A 1141B 1158 1159 
DIT VHD LI DIT 
UP LP VHD VHD 
\Vfole ° Mole Vole Vlole 
vercent percent percent percent 
OR. 62 OR 43 OR. #1 Os. 63 
03 3 03 02 
1.32 “Lig 33 1.31 
41 i) 19 , 4f 
¥i ‘ S4 
02 02 02 


02 02 . 02 


100.0 100.0 
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Runs 1141A and 1141 B were made on the same day. Runs 1) 

and 1159 were made 1 week later. The samples were withdray) 
from a standard 1-quart ICC container nearly filled with liquid, 
The technic employed was to set the cylinder-in an upright positioy 
and remove 0.4 mulliliter of liquid from the bottom valve by meay 
of a special liquid sampler, made from two small needle valves cop. 
nected by a short piece of 2.0 millimeter capillary tubing. . The enti 
quantity of liquid in the capillary was expanded into an evacuate 
glass cylinder, from:which one or more vapor samples were witb. 
drawn for analyses. Each analysis-refers to a separate liquid sam. 
ple computed on an entirely independent. basis. 
# The reproducibility obtained in an impure sample of butadiene js 
illustrated in table 6.. The extent to which the deviations given ar 
due to sampling and to the instrument. is difficult to tell. Erros 
from both sources should increase with the impurity. 


TABLE 6.—Crude butadiene 


Gom ponent Mole percent 





1,3-Butadiene........... 
1,2-Butadiene........... 
TTS ; 
1-Butene ey. 


n-Butane 
Methyl ethyl ether-.- 











® Limits represent the deviation from the average of duplicate analyses. 


This particular analysis is interesting in that it illustrates th 
procedure that must be followed when peaks are encountered for 
which no pattern is available. The sample was submitted -becaus 
it contained an impurity that interfered with both the infrared ané 
ultraviolet absorbtion analyses of butadiene. In addition to th 
usual impurities the pattern showed a set of peaks with primar 
contributions at 60, 45, and 31. These-masses can only result fron 
the fragmentation of n- or isopropyl alcohol, acetic acid, methiy! 
formate or methyl ethyl ether. Each of these compounds was ob- 
tained and its pattern determined; methyl ethyl ether. alone agri 
with the sample pattern. 

In the event that the pattern for the possible contaminants coult 
not have been obtained, the one actually present could have bee 
identified and its concentration estimated with a reasonable degr 
of accuracy. The ratio of 61 and 46 isotope peaks. to the 60 and 4) 
peaks showed only one oxygen atom to be present. This eliminated 
the formate and acetic acid. The two alcohols were improbable 0 
the basis of boiling points. The sensitivity of methyl ethyl ether 
the only remaining possibility, could be estimated from analogous 
compounds with very little error in the final calculations. 

The reproducibility and range of analyses in complex gas mixtures 

are illustrated in tables 7, 8,.and 9. 
The analyses given are representative of those obtained in copolyme! 
plants for aoe butadiene and for the noncondensableé vapors from 
the styrene’condenser, whereas table 9 is- representative of products 
of combustion of certain types of oil flames. 
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TABLE 7,—Recovered butadiene -° 


Component Mole percent 


1,3-Butadiene — F 70. 53 . 04 
1,2-Butadiene. - -- _ a 4.6 t+. 13 
Propadiene : . | 2. 4: 01 
Propylene.... . . ee { 10 
1-Butene___...... . 3.§ a 


2-Butene. _. . ‘ 215 +.09 
n-Butane ben : . ; 4 +. 
Pentadienes , eeestel : 01 
Pentenes_...- ‘ an ol . 05 . 001 
Eee ; | 5 . 02 


Vinyl cyclohexéne ° Ol 
Air . ; 01 
Carbon dioxide : jo +. 02 


Total 
iits represent the deviation from the average of duplicate analyses. 


TABLE 8.—Noncondensible gases from styrene condenser 


| 


Component Mole percent 


| 
1,3-Butadiene. - pe | 22. 96e+0.04 
Ethylene.....-- 73 08 
Propylene_-.... r ve 08 02 
1-Butene 
2-Butene. . 


n-Butane 
Pentadienes 
Pentenes._. 
Styrene 

Ethy] benzene 


PH PE 


Vinyl cyclohexene 
Nitrogen 

Oxygen 

Argon... 

Carbon dioxide 


Total 
represent the deviation from the average of duplicate analyses. 


TABLE 9.—Oil flame fumes 





Component Mole percent | 


.05 .Ol 
Ethylene......... -45 . 04 
Propylene ‘ 06. + .O1 
Acetaldehyde rae | .18 .02 
| 
Butenes aebehin | .02 +.01 
Benzene : Janeen .07 + .01 
i ndinccnmeedws splineenons .010 . 005 
. 06 .O1 
» | Ll 


Acetylene 


a " : --| 0. 48* +0. 05 


Carbon monoxide, ... 2. 91 . 03 
Oxygen...-. 3. 53 06 
Argon...... 10 .01 
nk nnccnimmancsboniiiincerenmcnnal 38 .02 


* Limits represent the deviation from the average of duplicate analyses, 
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Two magnet settings were required in taking the record for eae 
mixture. A magnet current of 0.5 ampere was used to sweep t! 
mass region frem 12 to 90 and a current of 0.7 ampere for the ‘regi 
from 90 to 150. 

The reproducibility of the instrument, as well as that of the sampling 
technic desertbed under table 5, is shown in table 10. 


TABLE 10.—Hydrocarbon- mixture 


Liquid sam ple /21/45 IT (3/22 
V apor-sampl 


component \fole percent Mole percent | Mole percent 
Proy ne | } a+) Ot 28 a+). 02 a+). 01 
}- Butadiene i a ] ». 5S : ! , + 01 
tal Bute ‘ Mi ¢ ; 4 03 
n-Butenes ‘ 7 1 5 4 D , 
lsobutene 
otal Butan 
n-Butar 
Isobutane 
otal Pentanes 
n-Pentane 
lsopentant 


‘ 


*® Limits represent deviation from the.average of 2 independent analyses of the san 


from @ vapor ré ervoir 

Kach double column headed by roman. numerals refers to a liqu 
sample removed from the container and completely evaperated into 
vapor reservoir... From this container, mass spectrometer inlet res 
vous designated as A and B were twice filled simultaneously, as gis 
by arabic numerals, k and 2. The limits designated by “‘a” are for | 
deviations from the average of the analyses of the vapor in A and in| 

The deviations from the average in each single column represent +! 
reproducibility of the instrument alone.. The deviations between 
and 2 represent reproducibility plus errors'in sampling the vap 
reservoir. Tho over-all spread under the roman numerals includes 
addition the lack of uniformity in liquid sampling. 

The impurities in tank helium from three. Government own 
plants are shown in table 11. 

TABLE 11.—Helium from various helium plants 


OTIS, KAN 
Impuritie 


Hydregen 
Nitrogen 


Hydrogen 
Nitrogen 


Potal 
NAVAJO, N 
Hydrogen 
Nitrogen 


Potal 
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The presence of hydrogen. was completely unsuspected; it was 
ater confirmed, however, by chemieal tests. .It was thought at first 
that thie iydrogen might have come from the steel storage cylinders, 
but helium from glass cylinders contained the same percentage of 
impurities. An analysis of the natural gas showed a trace of hydrogen; 
apparently it was concentrated along-with the helium in the purifica- 
tion process. ; 

The composition of the natural gas from: which 
obtained is shown in table 12. 


Texas helium is 


'T 


Taste 12.—Natural gas 


Mole 


Component percent 


Hydrogen 
Helium 
Neon 
Nitrogen 


race ® 
2. 24 +%).02 
123 Ol 


23.600 4 O5 


{rgon... 18+ .01 
COs, . 52 3 05 
Oxygen ‘ t 005 
Methane i. 16 4 02 
Ethane : 3.614 .02 
Propane | 2172 U2 
n-Butane bt .01 
Pentanes 36 4 02 
Hexanes 7+ U2 
Heptanes and higher. 


Total 


* Maximum less than 0.05 mole percent. 
> Limits represent deviation from average 
pendent an 


ilyses 
It is interesting to note that the neon to nitrogen ratio is appreciably 
higher than that found in the atmosphere, whereas the argon to nitro- 
gen ratio is slightly lower. 
The analysis of a special gas test mixture 
independent analyses were made by two computers. 


is shown ‘in table 13. 
The 
sample was also analyzed by the Gas Chemistry Section of the 


ie 
bureau. 


Four 


TABLE 13.—Gas mizture analyses 
10le percent 


FP) 588 


. 9.41 +0. 05 
wv - 02 
Unsaturated 
0.0 +0.02 


. 076 +. 002 
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Two magnet settings were used in taking the record, 0.15 amper 
below mass 12 and 0.5 ampere above. 

The separation of Nz and CO is notsharp. This is because the bag 
peak for both gases is 28. Separation must be made, therefore, by 
using the 14 and 12 peaks. For a 28 peak height of -100.00 the 
nitrogen 14 and carbon monoxide 12 peaks are 10.28 and 8.50, respec. 
tively; these adverse ratios result in a ten-fold magnification of the 
normal uncertainty. 


VI. TYPE'OF MATERIALS THAT CAN BE ANALYZED 


A sample to be analyzed successfully in the mass spectrometer 
should have an appreciable vapor pressure at room temperature 
Materials having vapor pressures as low as 10~* millimeter have been 
run in the instrument, but the peaks are too low for: precise calculs- 
tions. (In general, it is preferable that the vapor pressure exceed | 
millimeter for direct analyses.) A method will be described in suc. 
ceeding articles by which the analyses of rubbers, plastics, and mate. 
rials having no appreciable vapor pressure can be performed. 

The chemical composition of the material to be-analyzed is of little 
consequence. - For the sake of longevity of the filament it is desirable 
to keep water vapor low although samples saturated with water vapor 
at atmospheric pressure are regularly tested. Gas mixtures contain- 
ing small amounts of PH;, HS, SQ., and- chlorinated and fluorinated 
hydrocarbons offer no difficulty. 

The most difficult materials to separate are the stereoisomers 
The cracking patterns for cis- and trans-butene-2 are given in table | 
The fragmentation ratios are very close, the greatest difference being 
for mass'29. ‘The 58:29 mass ratios are sufficient to give a 1-percent 
separation in a mixture of pure butene-2. In a hydrocarbon mixture 
however, a separation of better than 2 percent cannot be expected, 
even under favorable conditions. The cracking patterns for all the 
pentenes have not been worked aut... 

The range of masses that can be resolved in the instrument covers 
all masses up to 150 for the slit settings ordinarily used. _ The isotopes 
of mercury can be distinguished but resolution. is not complete. In 
the analyses of hydrocarbons resolution above 150 is not essential 
as the sensitivity can be determined from any lower peak of sufficient 
size, and the pattern can be taken from any series of carbon groups 
on the record. The analyses of heavy hydrocarbons is limited more 
by vapor pressure than by mass. 


VII. ACCURACY OF ANALYSES 


A survey of some.1,700 determinations that have been run up to 
the present time leads to the conclusion that the precision of the mass 
spectrometer is high. Duplicate determinations in zeneral check to 
within a few tenths to a few thousandths of 1 percent, the usual devi- 
ation being a few hundredths of 1 percent. 

The absolute accuracy, in contrast, is not easy to determine. 4 
large number of determinations have been compared with the freez 
ape method. Agreements to within the limits of reproduci 
bility have been. obtained in every instance except oné, in which case 


the container from which the sample was withdrawn was nearly 
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empty: These results are not to be generalized, as freezing-point 
measurements can be made only on nearly pure materials. 

A number of checks have been madeagainst synthetic mixtures 
and with mixtures analyzed. by chemical methods, and agreements 
to within a few tenths of 1 percent have been obtained. Such results, 
however, cannot be used to ‘define the accuracy, as in no case was the. 
composition known to be within the limits of deviation. It is obvious 
that the limits of accuracy should not be defined by occasional extreme 
cases of agreement or disagreement. 

Errors in analysis enter from three different sources: (1) the in- 
strument, (2) the sample, and (3) the computations. 

The problem of obtaining a representative sample is far from siinple. 
The principal source of instrumental error results from a time-to- 
time variation in the sensitivity, or the cracking pattern or both. 
The effect of these variations can be largely overcome by checking . 
the pattern and sensitivity for the most abundant constituents in the 
mixture being tested and then correcting the patterns and sensitivitics 
of the other constituents. proportionally.. Thus, in the: analysis of . 
butadiene, the pattern error can be rendered negligible by running 
pure butadiene every day and correcting the patterns of-all the other . 
constituents therefrom. It should be pointed out that a daily check 
on the pattern is not imperative except in obtaining the precision re- 
ported in this paper. 

Instrumental errors may also enter through fractionation in the 
gas-inlet leak. When. capillary leaks are used, fractionation errors 
may be large; in the special leak used on the instrument this source 
of error is reduced to negligible proportions. 

Errors in computation, may enter through the use of incorrect 
patterns and sensitivities for the various constituents in the sample 
and through an improper interpretation of the pattern for the mixture. 

For the most part, instrumental errors are manifest by lack. of re- 
producibility between duplicate analyses, whereas computational, 
errors are indicated by the magnitude of the residuals, such as those 
shown in the last column of table 3 

In view of the number and nature of the various -contributing 
factors involved, itis impossible to define the accuracy of ‘a result 
obtained with the mass spectrometer by any given set of limits. In 
the analyses given in this ‘report there is sufficient reason to believe 
o the absolute accuracy is very close to the limits of reproduci- 

bility. Analyses have beer’ made, however, that contained mistakes 

judgment in interpreting the mixture pattern, and in: which in- 
correct calibrating patterns were used. In thése ‘instances ‘the re-- 
siduals were usually large. 
_In general, it can be stated that when the residuals are small the 
limits of accuracy are of the same order of magnitude as the limitsof 
reproducibility, and when the residuals are large the limits of accuracy 
may be correspondingly large. 


The authors are indebted to-H. W. Bond for preparing many of the 
hydrocarbons used-as standards, to Dorothy Thompson and Laura 
Edmundson for making many of thé computations shown, and to O. L. 
Parham for. operating the mass spectrometer. 

WASHINGTON, June 12, 1945 
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HEATS OF COMBUSTION AND ISOMERIZATION OF THE 
EIGHT C,H, ALKYLBENZENES 


By Walter H. Johnson, Edward J. Prosen, and Frederick D. Rossini 


ABSTRACT 


The heats of isomerization of the eight C,H,» alkylbenzenes were determined 
by measurement of the ratios of the heats of combustion in the liquid state of 
purified samples of these compounds by the procedure previously described for 
the hexanes, heptanes, and octanes. The data yield the following values for 
the heat of isomerization in the liquid state at 25° C, AH°, of'n-propylbenzene 
into each of the C,H; alkylbenzenes, in kilocalories per mole: n-propylbenzene, 
0.00; isopropylbenzene, —0.67 +0.16; 1-methyl-2-ethylbenzene, —1.93 +0.16; 
|-methy]-3-ethylbenzene, —2.48 + 0.19; 1-methyl-4- -ethylbenzene, —2.74 +0.26; 
1,2,3-trimethylbenzene, —4.83 + 0.22; 1,2,4-trimethylbenzene, —5:61 +0.16; 
|,3,5-trimethylbenzene, —6.00 +0.26. These data were combined with the 
value previously reported for the heat of combustion of n-propylbenzene to 
‘btain values for the heats of combustion of each of the C,H, alkylbenzenes in 
the liquid state at 25° C. 
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I. INTRODUCTION 


In continuation of the program of determining the heats of com- 
bustion, formation, and isomerization of hydrocarbons. of various 
types [1, 2, 3, 4],' calorimetric measurements have been made that 
yield values for the differences in the heats of combustion, or the heats 
of isomerization, of the eight C,H,, alkylbenzenes in the liquid state 
at 25°C. These heats of isomerization are combined with the values 
previously reported for the heat of combustion of n-propylbenzene 

| to obtain values for the heats of combustion of each of the C He . 


» beet in the liquid state at 25° C. 
—— 


' Figures in brackets indicate the literature references at the end of this paper. 
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II: UNIT OF ENERGY, MOLECULAR WEIGHTS, 
UNCERTAINTIES 


The unit of energy upon which the values reported in this papy 
are based is the international joule, derived from mean solar seconds 
and from the units of international ohms and international volts in terms 
of which certification of standards of resistance and electromotiy, 
force is made by the National Bureau of Standards. For conversioy 
to the conventional thermochemical calorie, the following relatioy 
is used [5]: 

1 calorie=4.1833 international joules. 


The molecular weight of carbon dioxide, the mass of which was used 
to-determine the amount of reaction, was taken as 44.010, from tly 
1941 table of atomic weights [6]. 

The uncertainties assigned to the various quantities dealt with i: 
this paper were derived, where: possible, by a method previously 
described [7]. 

Definitions of the symbols used are given in the previous papers 
[1, 2, 3, 4]. 

III. METHOD AND APPARATUS 


The same method and apparatus were used in the present investig- 
tion as in the work on the hexanes [1], heptanes [2], and octanes {4 
with the exception that a chronograph was tised to record the tin 
during the first 5 minutes of the reaction period. The chronograp! 
had two pens, one recording seconds and the other the taps of a tele. 
graph key. The telegraph key was tapped when the resistance of th 
platinum resistance thermometer in the calorimeter passed. certair 


previously selected values. 


IV. MATERIALS 


The compounds used in the present investigation were sample 
from the API-NBS series of highly purified hydrocarbons, which an 
being prepared through a cooperative undertaking of the America 
Petroleum Institute and the National Bureau of Standards.. A con- 
plete description of the purification and determination of the purity 
of these compounds will appear in other. reports [8, 9]. . The amounts 
of impurity in these samples, as détermined from measurements of 
freezing points made by A. R. Glasgow, Jr., E. T. Murphy, and A. J 
Streiff, at this Bureau; were as follows: n-Propylbenzene, 0.002 
+0.0008;: isopropylbenzene, 0.0005 +0.0002; 1-methyl-2-ethylbe: 
zene, 0.0024 +0.0007; 1-methyl-3-ethylbenzene, 0.0022 -+0.0015 
1-methyl-4-ethylbenzene, 0.0005 +0.0002; 1,2,3-trimethylbenzer 
0.0008 +0.0004; 1,2;4-trimethylbenzene, 0.0036 +0.0020; 1,3) 
trimethylbenzene, 0.0007 + 0.0003, mole fraction [8, 9]. 

As the manner of purification of these compounds [8, 9] was su 
as to leave substantially only close-boiling isomeric impurities in t 
respective compounds, it is calculated that in the worst case 1! 
measured heat of combustion would be affected by less than 0.00. 
percent because of impurities in the compound. 


V. RESULTS 


The experimental results of the present investigation are sumn 
-_ 


rized in table 1, which gives foreach of the eight comipounds the fo 
ing data: The number of experiments performed; the minimun 
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maximum values of the mass of carbon dioxide formed in the « 
bustion and of the calorimetric quantities /, K, U, AR,, Ar, 

the mean value of 2, and its standard deviation, in ohms per gr; 
of carbon dioxidé formed, as defined by equation 4 of referenc: 
The symbols in table 1 have the same significance asin the previ 
report [1] and the references there cited. 

In table 2 are given, for the eight C,H, alkylbénzenes, values 
the following: The constant B in ohms per gram of carbor dioxid 
given in table 1; B , which is B corrected to the ideal bomb react 
by the method of Washburn [10];? B?/B°, which is equal tothe ratio 
the heat evolved, per mole of hydrocarbon, in the ideal bo 
process at. 28° C for each isomer to that of n-propylbenze 
(+AU bk)» Al’,): the difference, between n-propylbenzene 
each isomer, in the heat of. combustion in the ideal bomb process 
28°. C; H; (liquid)—H, (liquid), the heat of isomerization. of 
propylbenzene into each isomer, at. 25° C and I atmosphere, for ¢ 
liquid state; and AHc?, the decrement in heat content accompa! 
ing the reaction of combustion of the hydrocarbon in the liquid sta 
in oxygen to form gaseous carbon dioxide. and liquid water, with 
the reactants‘and products in their the rmody namic standard states 
25°.C.. The value of the heat.of isomerization was obtained by mea 
of the relation -[4]: 


AUR) AU? AU R)n (1—B;/B,) 


For this calculation, the value of Al’,), at 28° C was take 
5209.2 int. kj/mole [3]. The value for the heat.of combustion « 
] 
given isomer was obtained by appropriately combining the heat 
isomerization with the heat of combustion .of n-propylbenzen 
given by the equations: 
C,H, (n-propylbenzene, liquid)+120, (gas)=9 CQO, (gas)+6H, 
(liquid). 
Ai c.16 5217.37 + 0.68 int. kj/mole=— 1247.19+ 0.16 keal/mo) 
The experimental data of the present and the previous investig: 
tions on thé alkylbenzenes [3] are being compared with the vr ta 
earlier investigations in a report [11] that will include selected“ 
values for the heats of combustion and formatien, for both the ic iq 
and gaseous states, at 25° C, of all the isomeric alkylbenzenes througi 
C,,H,, with general formulas for calculating values for the highe 
alkylbenzenes. 
? The Washburn correction is the same for all these mpounds as they are isomers, but accot 


the variation of the Washburn correction with amount ef sample burned As used here, the 
correction was modified to apply to 28° C and to the gase t zero pressure (instead of Latmosphe 
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THIN-WALLED ALUMINUM BETA-RAY TUBE COUNTERS 
By Burrell W. Brown and L. F. Curtiss 


ABSTRACT 


A description is given of aluminum-tube Geiger-Miller counters having a wall 
thickness of 0.004 inch and constructed from commercially available tooth-paste 
tubes. These tubes as furnished are in the hard temper and have uniform wall 
thickness, which enable them to stand complete evacuation without a tendency 

)collapse. By cepperplating the. aluminum it is possible to soft-solder fittings 

the tube to construct a counter that can be permanently sealed after filling and 
that will maintain its characteristics over a long period. The procedure reduces 
the cost of making this type of counter, as the aluminum tubes are inexpensive, and 
very simple operations are required in the construction of the counters. 


CONTENTS 


Ill. Copperplating the aluminum tube 
IV. Construction of the Geiger-Miller counter 


I. INTRODUCTION 


An economical. method of making thin-walled aluminum Geiger- 
Miller counter tubes for the measurement of beta rays would greatly 
increase the usefulness of this type of instrument in the field of radio- 
activity. Present methods require, the boring of metal tubes or rods ! 

}a thin wall. - This is a delicate.and expensive operation, and the 
thinnest walls obtainable are approximately'0.007 inch. Thin-walled 
counters have also been made by drawing glass tubing to a thin-walled 
section and silvering the inside.? This procedure does not give a very 

iform wall, which is usually of a thickness greater than that corre- 
sponding to 0.007-inch aluminum. Other methods require the use of 
wax seals, which do not yield tubes of very long life as the organic 
vapors commonly used ‘in filling counters attack most waxes. 


II. THIN-WALLED ALUMINUM TUBES 


It has recently been noted that thin-walled aluminum tubes quite 
suitable for the construction of counters are commercially available. 
4 typical tube is % inch in diameter and 4% inches long, with walls 
0.0044+ 0.0004 inch in thickness. The walls are uniform throughout 
: J R. Dunning and 8. M. Skinner, Rev. Sci. Instr. 6, 243 (1935); J. Halpern and O. C. Simpson, 8, 172 
vA H. Barnes, Rev. Sci. Instr. 7, 107 (1936); J. C. Wang, J. F. Marvin, and K. W. Stenstrom 13, 81 (1942). 
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the length and circumference within the limits indicated. These-{ 
are made in quantity as containers for such materials as tooth+ 
and are furnished with one end fitted with the conventional screy 
and the. opposite end open to the full diameter of the tube. Mad 
mass-production methods these tubes are inexpensive and of 

uniform dimensions, and therefore are suitable as a component , 
beta-ray counter. They are in the hard temper, which enables | 
to withstand complete evacuation with no tendency to collaps 


yas 


III. COPPERPLATING THE ALUMINUM TUBE 


It is of advantage to construct counters that can be filled 
sealed permanently so that they will maintain their charact: 
over considerable periods of time. This requires the eliminatio 
wax seals, gaskets, and similar arrangements. To accomplish 
with the aluminum-tube counters, it has been found desirable to p 
the aluminum with copper to permit the use of soft solder in jo 
the parts. There are several methods by which an adherent plat 
of copper can be applied to aluminum. For the benefit of thos: 
may desire to construct such counters, the method of: plating fi 
most convenient is described. 

The first step is to clean the aluminum surfaces to be plated 
strong solution of potasstum hydroxide, immersing the alun 
until a smooth film of water will spread over it uniformly, indicat 
that the aluminum surface is clean. ‘The next step is to anodiz 
aluminum surface. After rinsing all the potassium hydroxide solut 
from, the aluminum, it is used as an electrode in a solution of 10- té 
percent oxalic acid to produce an anodic-film on the surface. | 
convenient to anodize two pieces together by using them as elect: 
connected to the alternating-current supply through a variable-vol 
transformer. Initially; the voltage applied should be of the ord 
| or 2 volts, and this may be increased to around 50 volts as the a 
film develops. ‘Care should be taken to hold the voltage low enouw 
at all times-to avoid heating of the solution by the current. — A! 
15 to 20 minutes the aluminum. is removed, rinsed in cold water 
etched for a few seconds in a 1-percent hydrofluoric acid solution 
a slight evolution of gas is apparent. It is then immediately rins 
cold water, and the surface is ready for plating. The. usual coy 
sulfate solution, to which a small amount of sulfuric acid has } 
added, is used fer plating. The. plating of a layer of copper a 
thousandths of an inch in thickness requires about 20 minutes 
constructing the counters, both ends of the tubes are plated \ 
copper. If copper is deposited over a larger area, than is desired 
excess may be removed at any time by etching with strong nitric 
which has little effect on the aluminum. 

The etching of the anodic coating. with hydrofluoric acid may 
eliminated. if care is taken to develop a thin coating. in the anodiz 
process. ‘The proper thickness is obtained by interrupting the an 
ing process when the potential reaches about 2 volts. Some diffi 
may be encountered in determining the proper conditions to secu! 
adherent coating of copper by use of this alternative method. 
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FIGURE 3 Completed counter 
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Ly. CONSTRUCTION OF THE GEIGER-MULLER COUNTER 


assembly of the counter is shown .in cross section in figure ‘1. 
circular end of the tooth-paste tube, Bol is closed by the copper 
ishing, B, which is soft soldered -to the plated aluminum tube, as 
eated. The bushing is fitted with a \4-inch Kovar tube, A, silver- 
lered in ar axial hole, as shown. A short glass tube sealed to K 
ovides a support for one end of the central wiré, W, and also carries 
side tube for evacuating, filling, and sealing the counter: The glass 
Rube Is continued inside of AK to form a sleeve around the central wire.’ 
‘he opposite end of the aluminum tube is closed by.a copper plug, P, 
ered into the small tubular opening of the tube. The plug is 

| te take a small piece of Kovar wire, which carries a glass bead 
erve as a support for the ceritral wire at, this end. The plug 

are soft soldered in one operation to the aluminum tube, 

viously copperplated at this point. 

‘he order of assembly found convenient is first to silver solder tube 
nto copper bushing, B. The glass tubing is then sealed into A and 
bushing, B, soft soldered into tube, 7. The central wire system is 
nade up and lawered into 7 to protrude through the screw-cap 
Copper plug, P, is soft soldered to the short Kovar wire and then 
The Kovar wire,., is then pulled into position at the opposite 








Cross section of assemi ly 


Land sealed into the glass. The dimensions are such that there is a 
ate tension on the-spring, S. If.the dimensions are carefully 
sen this spring may be eliminated. 
he performance of the finished counter depends, of course, on the 
of filling used. In the counter described here, the usual alcohol- 
m mixture has been found -to yield plateaus and thresholds similar 
other types of Geiger-Miiller counters containing this mixture. 
slope of 0.1 percent per volt is readily. obtained. <A curve showing 
iation of counting rate with voltage applied to the counter is 


1in figure 2.:° This curve is maintained for several months, which’ 


ng as. the counters were under observation. [Experience indi- 
tes that a- counter will leak occasionally as a result of imperfect 
ering. -This can be quickly repaired. In a few cases it was also 
nd that the aluminum tooth-paste tubes have a small leak, probably. 
»pin holes. Perhaps 5 perceint of the tubes may have this defect, 

ch can be detected by preliminary evacuation. 
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Ficure 2.—Varialion of counting rate with voltage applied to counter. 


A completed counter, ready for filling, and an empty tube a 
shown in figure 3.° Aluminum tubes’ of other diameters but 
approximately the same wall thickness and that do not collapse « 
evacuation are also available commercially. 


WasuinecTon, May 4, 1945. 
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METHOD OF cance CURVILINEAR, 
INTERPOLATION: 


By, William"J. Taylor? 


ABSTRACT 


This report describes a simplified method of computing Lagrangian coefficients | 
for curvilinear interpolation, which may be used when tables of Lagrangian 
efficients are not available or when tables are available but the coefficients are 
ot tabulated for the exact fraction of the interval to which the interpolation is 
0 be made. 
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I. INTRODUCTION 


The common formulas of curvilinear interpolation,’ such as those 
pf Newton, Gauss, Stirling, Bessel, and Lagrange * when carried to 
ompletion, are all equivalent in that they approximate the function 
y a polynomial of suitable degree (the Ble of points required is 
ne greater than the degree of the polynomial assumed). All of these 
ormulas except that of Lagrange require, as a first step, the compu- 
ation of a ‘table of differences” for the function, to an order. equal 
0 the degree of the polynomial. Lagrangian interpolation, on the 
ther hand, is carried out by multiplying eachi of the tabulated values 
bn which the interpolation is to be based by a suitable coefficient— 
he Lagrangian coeflicients—and taking the sum of these products. 
his “cumulative”? multiplication: can be carried out as a continuous 
peration on a modern.c diction machine, and when such a machine 

available Lagrangian interpolation is usually more rapid than the 
ther methods referred to. 

Very extensive tables of Lagrangian lntermeiation coefficients are 
low available * (references to other, less‘ co omplete, tables are given in 
us work). 

investigation was performed at the National Bureau of Standards as part of the work of the Ameri- 
leum Ir stitute Research .Project 44 on the ‘Collection and Analysis of Data on the Properties of 
re h Assoc iate on the American Petroleum Institute Research Project 44 at the National Bureau 
cussion in this re ort is restricted to interpolation in tables in which a function is tabulated at 
paced intervals of the independent variable. 
cas tt bo a eens Tracts No. 2, Interpolation and nymerical integration (G. 


Ables of Lagrangian interpolation coefficients, Mathematical Tables Project, Works Project Admin- 
BAe leral Works Agency. (National Bureau of Standards; Columbia Unive rsity Press, New York, 
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In this.report there is described a simplified method of comput 
the Lagrangian coefficients, which may be used when tables of | 
vrangian coefficients are: not available, or when tables are ay 
able, but the met or are not tabulated for the exact fraction 
the interval to which the interpolation is to be made. ‘This met 
requires only a single division for the computation of each coeflix 
for interpolation on the basis of any number of tabulated valu 
The coefficients, as thus computed, are not normalized—they may 
normalized in the usual way, although it is usually more conver 
to interpolate by means of the non-normalized Lagrangian coeflici 

Lagrangian interpolation, with the use of tables of Lagrar 
coefficients (see footnote 5), or by the metliod described in this rep 
is recommended as a rapid and generally useful method of inter 
tion in-those tables in which a function, or property, is tabulat 
several equally: spaced values of the independent variable. ' 
number of tabulated values, or the degree of the polynomial 
which the interpolation is based must be chosen to represent 
function with suflicient-accuracy in the interval in question. 


II. GENERAL METHOD 


The form. of the modified Lagrangian interpolation formula \ 
stated, and then the derivation of the formula will be given. 

Let the (n+1) tabulated values of the function y(z) on whicl 
interpolation is to be based be wy, w. .—s let 
corresponding values of the independent, vi ariable equally spac 
inte rvi als of h be To, 21 Ms | ‘I he (n+1) values ( 
determine a polynomial of the nth degree in z, P(x), which is an apy 
imation. to the function y(z) in the range z, to z,. The value of P 
and therefore, approximately, of y, at a given value of.z, is give! 
the modified Lagrangian formula 


where the coefficients, ay, 


a, 


(x7,—z)/h} is the absolute value r,.—z)/h, and (7) 1s the 


coefficient 


The signs of the’ coefficients, @,, rmined as follows 
hes between z,, and z, 
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coefhicients, ca 0) Am, Bmiity Bm+s) - , are positive, 
tS, @ » Am-1, Aur2, Omie; are negative. 


rus is thus of the type 


e] 


t 


r.Or INTERPOLATION 


are proportional to the Lagrangian interpolation 
In fact, the sum of the a,’s is 


8.-@ 
not normalized. 


‘SSiON 


a check on the calculation of the Cn 


in erpolations to the same fraction of an interval are to 


mio converment to put eq 1 in the form 


P(x) A Uky 


h ( 


are the normalized Lagrangian coefficients, 


have been tabulated clsewhere see footnote ). As usecl 
indicates that the algebraic sum of the coefficients 


“normalized”’ 


be derived in the following way: P(2) may be 


esults may 





by definition) and the expression on the right of eq 7 are both 


ls of the nth degree in z, and are equal to y, when 2 
But two polynomials of, the ath degree which 
points must be-identical, which establishes eq 7. 


z,, for 


fi 
at ~t- | 


o! eq 7 vields 
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r 

The expression in brackets under the summation sign is equal y 
absolute value to the binomial coefficient, (7). Therefore if the sigy 
of the terms are temporarily disregarded, as indicated by the + sin, 
eq 8 may be rewritten as 


Pa) =| = it |@—2y a} |S] £(2)/ (ex—2)/h| | 


Inspection of eq 7 shows that the signs of the terms follows the nj 
already given for the signs of the coefficients a;. In view of the defi 
tion of ay, eq 2, therefore, eq 9 becomes 


1 n | P | n 
a) =| TI |(a—z,)/h [2a 
TU k=O | Jk=0 


In the special case in which y,=1, fork=0,1, . . . ,n, the polynomi 
P(x) must have the value unity for all values of z. It follows that 


E II (x—2,)/h >a, 5 
Ti k=O | k=0 


Equations 10 and 11 lead immediately to eq 1, and eq 11 is equiv: 
to eq 4. 
III. NUMERICAL EXAMPLE 


Given the following (equally spaced) values of the descend 
exponential function, e~*, find by interpolation the value of ¢* «& 
z=0:54316: 


y=c* 


Yo= 0.594 520 548 0 
Y= 0. 588 604 969 7 ‘ 
Yy2= 0.682 748 252 4 
v3= 0.576 949 810 4 
¥s= 0.571 209 063 8 





Solution: The interval, A, of the independent variable, z, is | 
The values of the binomial coefficients for n=4 are 1, 4, 6, 4, and 
The value z=0. 54: 316 lies bétween 2.=0.54 and x3;=0.55.° Then! 
values of the a,’s, or the non-normalized Lagrangian coefficients, ! 
given by eq 2: 

ao= + FG) Go— r)/h|=+1/2.316=+ 0.431778 9 
ay - (f)/|(ay— 2) /h| = —4/1.316=— 3.039 513 7 
a,= + + Gilles x) /h| = +6/0.316= +18.987 341 8 
a3= + ($)/|(z3—2) /h| = +4/0.684=+ 5.847 953 2 
a,=- (iter z)/h|=—1/1.684=— 0.593 824 2 


4 
>) d= 21.633 736 0 


k=0 
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check may be obtained from eq 4: 


at/ it|(ex—2)/h| +24/(2.316) (1.316) (0.316) (0.684) (1.684) 
k= ; 


=21.633 736 0. 


The value of P(x) at r=0.54316 may now be calculated from eq 1 
he cumulative multiplication involved may be performed as a 
pntinuous operation on the calculating machine): 


r 4 
P(z)=>5 axys/ 33 ay 
k=0 k=0 


a 4 ers) 
~ (21.633 736 0 


=0.580 909 673. 


The true value of y=e-*, at z=0.54316, correct to 10 decimal 


paces, 18 


y=0.580 909 674.4. 


When the same coefficients are to be used several times, it is con- 
nient to compute the normalized Lagrangian coefficients, or the . 
3, from eq 6: 
A,= +0.019 958 592 
A,;=— .140 498 789 
A,;=-+-. .877 672 807 
A;=+ .270 316 380 
A,=—. .027 448 990 


4 
>> Ai= ‘1.000 000 000 


f=0 
hese coefficients differ by one or two units in the last place from the 
rrect-values (see p. 232 of the-reference in footnote 5) because the 


's weré computed to only seven decimal places. ‘The value of P(z) 
z=0.54316 may now be calculated from eq 5: 


4 
P(z) =>5 Ayr 
k=Q 


0.580 909 674: 


Wasninetron, March 3, 1945. 





